GAMOW VECTORS AND BOREL SUMMATION 



O. COSTIN, M. HUANG 

Abstract. We analyze the detailed time dependence of the wave function 
ijj{x,i) for one dimensional Hamiltonians H = —9^ + V{x) where V (for ex- 
ample modeling barriers or wells) and 5/1(0;, 0) are compactly supported. 

We show that the dispersive part of jp{x,t), its asymptotic series in powers 
of t~^/^, is Borel summable. The remainder, the difference between ip and the 
Borel sum, is a convergent expansion of the form ^'^—oOk^ki^)^ where 
Fj. are the Gamow vectors of H, and 7^. are the associated resonances; generi- 
cally, all are nonzero. For large k, 7^. ~ const ■ k log k + k'^-n'^i/i. The effect 
of the Gamow vectors is visible when time is not very large, and the decompo- 
sition defines rigorously resonances and Gamow vectors in a nonperturbative 
regime, in a physically relevant way. 

The decomposition allows for calculating t/i for moderate and large i, to 
any prescribed exponential accuracy, using optimal truncation of power series 
plus finitely many Gamow vectors contributions. 

The analytic structure of ip is perhaps surprising: in general (even in simple 
examples such as square wells), ip(x,t) turns out to be C°° in t but nowhere 
analytic on R+. In fact, i/; is t— analytic in a sector in the lower half plane and 
has the whole of a natural boundary. 

Extension to other types of potentials, for instance analytic at infinity, is 
briefly discussed, and in the process we study the singularity structure of the 
Green's function in a neighborhood of zero, in energy space. 
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1. Introduction 

Resonances play a major role in the physics of metastable states and their decay. 
From a mathematical standpoint, there is a good number of definitions of resonances 
and resonant states. In most approaches, they are based on the properties of the 
scattering matrix, on Gelfand triples (rigged Hilbert spaces), or on the complex 
analytic singular structure of the Green's function beyond the spectrum of the 
resolvent. The pole positions of the Green's function, "resonances" are pseudo- 
eigenvalues, and their residues (Gamow vectors) are pseudo-eigenvectors of the 
Hamiltonian with "purely growing" conditions at infinity. There is a vast literature 
on the subject, see e.g. the concise overview [13j and the references therein. See 
also [9 for a surprising consequence of resonant states, and for a clear description 
of the physical relevance of Gamow vectors. 

By and large, the different mathematical definitions provide equivalent objects. 
However, there are conceptual difficulties in all rigorous approaches, and these lie 
in connecting (a) the mathematical definition, (b) the natural properties of the 
underlying quantum Hamiltonian, and (c) the physical phenomenon. In fact. How- 
land's Razor, a principle so dubbed by B. Simon, cf. [_! 9J , states that no satisfac- 
tory definition of resonance can depend on the structure of a single operator on an 
abstract Hilbert space. Slightly oversimplifying the argument, the reason is that 
the analytic structure of the Green's function, or of quantities obtained through 
dilation-analyticity, needed in most approaches, are by no means unitarily invari- 
ant. Unitary invariance plays of course an important role in quantum mechanics 
since observables are represented by self-adjoint operators on Hilbert spaces, the 
isomorphisms of which are precisely the family of all unitary transformations. 

A concise and very illuminating critical analysis of the various mathematical 
attempts at rigorous definitions is found in [19j . 

1.1. Resonances and asymptotic expansions. We note however that many 
relevant physical quantities are not and need not be defined in a unitarily-equivalent 
way. As already mentioned, resonances are used in measuring the time decay of 
the probability distribution in physical space. In any interpretation of quantum 
mechanics, L'^(M'^) plays a distinguished role, when is a representation of the 
space where we, and macroscopic apparatuses, lie. 

A definition based on time behavior is natural to the underlying physics and 
avoids Howland's razor since it rests on (i) a particular representation of the 
Hamiltonian-as an operator on of our R^, (ii) on a second observable, say 1^, 
the characteristic function of the set ^ € and (iii) on a specific mathematical 
question-the time decay of This triad is not (at least not manifestly) 

a property of a single operator. Nonetheless, L^(]R^) and dependence on time are 
canonical objects in analyzing scattering or decay problems. 

At the present time however rigorous definitions based on time behavior only 
exist in a perturbative regime, [20], [14]; see also below. 

In this paper, for compactly supported potentials in one dimension, we show that 
the difference between the wave function and the Borel sum of its asymptotic series 
in powers of t^^^^ is a convergent expansion in Gamow vectors. The resonances 
thus defined turn out to be independent of the initial condition. Gamow vectors 
are not functions; neither is the Borel sum (see SjSj) of the power series. The 
expansion is valid uniformly on compact sets instead. 
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The representation as a Borcl summed series plus Gamow vectors expansion is 
valid not only for large t, but, in fact, simply for t > 0, though the convergence 
rate of the whole expansion is rapid enough only if t is not too small. 

After completing a manuscript we found that decompositions in energy space in 
terms of Gamow vectors and a continuous part have been proposed in the physics 
literature, see [8J, to our knowledge without completely rigorous, mathematical, 
proofs or study of Borel summability, and with a different interpretation and sug- 
gested physical meaning; cf. Note[5]below. Without Borel summability, uniqueness 
of a decomposition in terms of a continuum integral and a sum of exponentials 
generally does not hold, see ^JS) 

The power series expansion roughly corresponds to the decay of a free 

particle 0. Indeed, if time is very long and the point spectrum of H is empty, then, 
eventually, the overlap between the wave function and the support of the potential 
becomes negligible. The specifics of the potential are seen while the particle has a 
fair probability of it being near the potential. This is why it is natural to subtract 
out the power series, "free" decay. But, generally, the series has zero radius of 
convergence^. 

If parameters are such that a resonance (complex generalized eigenvalue, [14]) is 
at a small distance e to the spectrum of H, the setting is called perturbative and 
there is a time scale, roughly given by e~'^* 3> t^^^^ , during which in a finite spatial 
interval, the decay of the position probability follows an exponential law. This 
corresponds to a transient, metastable state. The Gamow vector corresponding to 
such a resonance describes the wave function on increasingly larger spatial regions, 
see [9], §9. Only metastable states with long enough survival time are captured 
however in this way. (Rigorously speaking, we are dealing with a double limit, 
in which time goes to infinity and an external parameter goes to zero in some 
correlated fashion.) Borel summation provides an exact representation for all i > 0, 
as well as practical ways to calculate the wave function for times of order one, see 
§5.1) the influence of resonances which are not necessarily close to the spectrum is 
measurable. 

For showing Borel summability, perhaps the most delicate part is the analysis 
of the Green's function in the fourth quadrant in the energy parameter, where 
infinitely many poles recede rapidly to infinity; sharp estimates are needed in order 
to control a needed Bromwich contour integral. 

Extension to other potentials with sufficient analyticity and decay is discussed 



2. Setting and main results 
We consider the one-dimensional Schrodinger equation 

d h d"^ 

ih—Tp{x, t) = t) + V{x)il){x, t) 

where: 



^The influence of the potential —however distant- is still present in the "initial state" at some 
very late time ti ^ 1, from which the almost free particle decays; the state at ti is responsible for 
the generic disappearance of the zero energy resonance. 

■^If the potential is unbounded, such as a dipole V{x) = Ex, then the power series may be 
identically zero, see |21j and references therein. Another exception is y = 0, for which the 

asymptotic t~^^^ series converges on compact sets in x. 
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(a) The nonzero potential V is independent of time, compactly supported and 
on its support. {V is allowed to be discontinuous at the endpoints provided that 

it is one-sided at the endpoints.) 

(b) The initial condition ipoix) is compactly supported and on its support. 
We normalize the equation to 

(1) i^^Hx,t) = -^^{x,t)+Vix)^ix,t) = {H^){x,t) 

where supp(T^) C [—1,1]. Under the assumptions above, we have the following 
results. 

Proposition 1. For large t, the wave function ^p{x,t) is 0{t^^^^) (in the generic 
case of absence of zero energy resonance (ID] . ip{x,t) — 0{t~^^'^)), and 'il}(x,t) has 
a Borel summahle asymptotic series ^{x^t) in powers of t^^/"^. 

We denote as usual by LB the Borel summation operator. Let t~^^^ip{x,t) = 
CB'>p{x,t) where ip{x, ■) is bounded. As seen below, 'ip{x,t) — t~^^^ip{x,t) is nonzero, 
and is a convergent combination of Gamow vectors, the residues at the poles of the 
analytic continuation of the resolvent of H. 

Let {Ek}k=i,...,N be the eigenvalues of H and at be the corresponding 

eigenf unctions. (We convene to set = if these two sets are empty.) Let 
also 7fc,Re7fc > be the generalized eigenvalues (resonances) corresponding to the 
Gamow vectors Tk{x). 

Theorem 1. (i) For all t > we have 

N oo 

(2) V(x, t) - t-^^'^{x, t)^Yl bkM^)e-''''' + J2 9krk{x)e-'^'' 

k=l k=l 

The infinite sum in 0) is uniformly convergent on compact sets in x -rapidly so if 
t is large. (The coefficients bk and gk depend on -0 and typically gk =/= for all k.) 

(ii) ipk{x),Tk{x), ip{x,t) are twice differ entiable in x. 

(Hi) We have 

(3) 7fc ~ const ■ fclogfc + fc^7r^i/4 as k +oo 

(Higher orders depend on V, see Proposition\^) The 7^ are independent of ipo, 
and the constant depends on the endpoint behavior of V . 

The series in ([2]), though valid for all t, converges poorly if t — > 0: this is not the 
regime it is intended for. 
Let 




where E„ is the n-exponential integral and argit e (— tt, 0) U (0, tt), [T] 



■^There seem to be inconsistent definitions in the literature. We use _1 ; since the definition is 
not spelled out in one place, we state it again: E„(2) = f^ t^"e~^*dt for Re^ > 0, analytically 
continued to C \ R~ , and extended by continuity to the two sides of the cut. 
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Proposition 2 (o(e ^^*) accuracy, for arbitrary M). For any M there exists (ex- 
plicit) m and mi < m, so that 

N mi rn 

(5) ij{x, t)=J2 &fcV'fc(a;)e~^'=" + J2 9kTk{x)e-^''* - ^ rkS{%t) + i^Mix, t) 

k=l fc=l k=l 

Here {'Jk}k<m the poles of the Green's function (resonances) on the first and 
second Riemann sheet with \^k\ < M , {7fc}fc<mi are ihe subset of them on the first 
Riemann sheet, r^ are the corresponding residues, and i/jm differs by o{e~^*) from 
its (relatively explicit) power series in l/t^/"^, optimally truncated (see ^5.1\) . 

Note 1. (i) We see that, as exponential contributions, only the resonances on the 
first Riemann sheet appear but both sheets contribute to the dispersive part. 

(ii) The expression ([4]) is not analytic on the Riemann surface of the log: it has a 
jump on K"*", compensated by an opposite jump of ^m- These jumps are mandated 
by least term summability requirements. 

Corollary 3. Any number of resonances can be calculated from tp{x,t), if tp is 
known with correspondingly high accuracy. Conversely, ■0 can be calculated in 
principle with arbitrary accuracy from the contribution of a finite number of bound 
states, resonances, exponential integrals and optimal truncation of power series. 

Note 2 (Analytic structure of ipix^t), in t e K+). It follows from the proof that 
In |(7/i;r/i;| = 0{,Jpk). Thus, since Lp is manifestly analytic for Re t > 0, it follows 
immediately from ^ that ^j: is C°° in t. Now, since Im — 7^ ~ —k^, near R+, 
?A equals a function analytic in the right half plane (the Laplace transform) plus a 
lacunary Dirichlet series, convergent for Im t < (the "heat-like" direction). For 

generic x, the coefficients of the Dirichlet series are bounded below by e~™"'''(^)V^ 
(all functions involved are of exponential order 1/2; the lower bounds follow rela- 
tively easily from the proofs, but we omit the details). Then the Dirichlet series 
does not converge past K"*"; general theorems on lacunary series, see e.g. |15| imply 
then that M"*" is a natural boundary. See also Proposition[8] Using similar estimates 
it can be checked that the Taylor coefficients of V' at a point to behave roughly like 

k 

, showing once more that there is no point of analyticity on M+. 



^eHl \n k 



This is another way to see the contribution of the Gamow vectors to the properties 
of ip. (More details about this are part of a future paper.) 

3. Proofs of Main Results 

3.1. Integral reformulation of the problem. H satisfies the assumptions of 
Theorem X.71, llSj v. 2 pp 290. Thus, for any t, tp{t, •) is in the domain of -d^/dx"^. 
This implies continuity in x of ip^tjx) and of its t— Laplace transform. It also 
follows that the unitary propagator U{t) is strongly differentiable in t. Existence of 
a strongly differentiable unitary propagator for ([T]) implies existence of the Laplace 
transform 

^{x,p)=J^ e-P'il>{x,t)dt^ (^j^ e'P*U{t)dt^Mx) 
for Rep > 0. Taking the Laplace transform of ([1]) we obtain 

(6) ip^{x,p) - i^poix) = --^ip{x,p) + V{x)%p{x,p) 
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where "00(2;) is the initial condition. Treating p as a parameter, we write ip{x,p) = 
y{x;p) —: y{x), and obtain 

(7) y"{x)^iV{x)~ip)yix)^tMx) 
where y{x) £ L^(R). The associated homogeneous equation is 

(8) y"ix)^{Vix)-ip)yix) 

If y+{x),y-{x) are two linearly independent solutions of ([8]) with the additional 
restrictions (and the usual branch of the log) 

y^{x) — e^^/^^ when x > 1 

(9) y_(a;) = e^'^^^ when x < -1 

then, for Rep > 0, the solution of ^ (or equivalently of (I?])) is 

(10) ij;{x,p) = (y^{x) y+{s)^jjQ{s)ds - y+{x) 2/_ (s)V'o(s)ds 

V J+oc J~oo 

where the Wronskian Wp — y-^-{x)y'_{x) — y^(x)y'^{x) is easily seen to be indepen- 
dent of a:. 

As we shall see, this solution is meromorphic in p except for a possible branch 
point at 0, and for fixed x it has sub-exponential bounds in the left half p-plane 
(when not close to poles). The function "0 is the inverse Laplace transform of "0, 
and it can be written in the form 0(a;,t) = lag-i^ ^{x,p)eP*dp. We show that 
the contour of integration can be pushed through the left half plane; collecting the 
contributions from poles and branch points, the decomposition follows. 

Note 3. The domain of interest in p is a sector on the Riemann surface of the 
square root, centered on R"*" and of opening slightly more than 2n, which, in the 
variable \f—vp translates into a sector of opening more that tt centered at \f^i. 

3.2. Analyticity of ^ on the Riemann surface of y^. We start with the 
analyticity properties of '(/'■ The more delicate analysis of the asymptotic behavior 
of the analytic continuation of on the Riemann surface of the log at zero is done 
in H3.4I The existence of a square root branch point at zero is typical in this type 
of problems. For our analysis, in proving Borel summability, we need to show that 
■0 is meromorphic in ^/p. 

Proposition 4. ip{x,p) is meromorphic in p on the Riemann surface of the square 
root at zero, Ci/2;o o,nd zero is a possible square root branch point. 

Proof. This follows from the following simple argument. Note first that continuity 
of y and y' imply the following matching conditions: 




Consider now the solutions /i and /2 of ([S]) with initial conditions /i(— 1) = 
1, /{(—I) = and /2(— 1) = 0, /2(— 1) = 1- By standard results on analytic 
parametric-dependence of solutions of differential equations (see, e.g. [H]), we see 



GAMOW VECTORS AND BOREL SUMMABILITY 



7 



that /i and /2 are defined on R and for fixed x they are entire in p. We note that, 
by construction, the Wronskian [/i,/2] is one. Then, 

y+{x) = Cihix) + C2/2(x), y.{x) = C^h{x) + C4/2(a;) 

where 

Ci = V^e-v^(/2(l)-/^(l)) 
C2 = -v^e-v^(/i(l)-/{(l)) 
C3 = -v/^e"^^ 



'ipe 



Furthermore, 



(11) W, = (*P(/2(1)) + /{(I)) - V^(/i(l) + /^(l))) 

Thus y± and Wp are analytic in Ci/2;o with a possible branch point at zero. The 
same follows for ip, by inspection, if we rewrite its expression as 



(12) ^{x,p)^ — {y^{x){ / y+is)Ms)ds+ e"^-'P'Ms)ds 



+ 00 

y+{x) ( / y-{s)ipo{s)ds + 



□ 



3.3. The poles for large p in the left half plane. To effectively calculate the 
asymptotic position of poles as p — > cx) in the left half plane, we need a more 
convenient choice for /i , /2 . In the previous subsection they were chosen to be 
analytic in p. Here we choose a new pair of /i , /2 for which the asymptotic behavior 
as p — !■ 00 is manifest. 

Note 4. It is straightforward to check that if fi{x) and f2{x) are solutions of 
([H]), and their Wronskian Wf-p = [/i,/2] is nonzero, then in the decomposition 
y+{x) = Cifi{x) + C2f2{x), y-{x) = C^fiix) + (74/2(2;) we have 



C2 = 

C3 = - 

Ca = J- 



-tp 



-IP- 



-ip- 



(/2(1)-/^(1)) 



(/l(l)-/((l)) 



-ip- 



ip- 



Wf;p 



(/2(-l)+/^(-l)) 

(/i(-l) + /{(-!)) 



Furthermore, Wp = [y+,y-] = (C1C4 - C2Cz)[fij2] is given by 

-2V=ip / 

(13) Wp = -— a(A(l)/2(-l) - /i(-l)/2(l)) - /{(l)/^(-l) + 



/((-1)/2(1) - \/^(-/{(-l)/2(l) - /{(l)/2(-l) + /l(l)/^(-l) + /l(-l)/^(l) 
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Proposition 5 (WKB solutions). In S+ = {p : Ke{y/—ip) > 0} there exist two 
linearly independent solutions of (0) of the form 



(14) 
(15) 



-fl--^ r VtV{s)ds + -9iix) 

V ^Vp Jq p 



where gi(x), g[{x), g2{x), g2{x) are bounded in p as p oo in 5+. A similar state- 
ment holds S- = {p : Iic{\/—ip) < 0} 

Note 5. This is in a sense standard WKB; however, since details about the regu- 
larity of the terms expansion are needed we provide a complete proof. 

Proof. We will only prove the conclusion for gi, since the proof for 172 follows 
analogously. Substituting into ^ , we obtain the equation for gi : 



-ip 



(16) g'l{x) - 2^/~Apg[{x) - V{x)gi{x) + i- 
We rewrite this equation as an integral equation for g'^: 

(17) g[{x)^e^^P^ 

[u)du- 



V{s)ds - V\x) = 









Vis) [ 




Jo 



ds 



where = 1 if — 7r/2 < argp < 37r/2, and xq — —1 if — 57r/2 < argp < — 7r/2. 
Note that |e~2V-jp('<-a:) | ^ \ fQj. g between xq and x. Using integration by parts 
we obtain 

(18) g[{x)^-^-^ (v{x) j\[(u)du-e-^^P^^'>-^W{x^) j^\[{u)du^ 
+ £ e-2/-P(--) (^^'(,) + V{s)g[[s)^ ds 

- ^ V{u)du - V'{x) - e-2v^(-o-^) (^j^" v{ti)du - y'(xo)^^ 

+ 1 r e--^^i-'-^) {V{s)-V"{s))ds 



For large p, under the norm ||/|| — sup^g[_;^^] \ f{x)\ the above integral equation 
is easily seen to be contractive inside the ball 



^ sup 

-l<a;<l 



V{x) 



V'{x) 



V"{x) 



Therefore g'i{x) and gi{x) — g[{u)du are both bounded in p as p 00. □ 



Remark 6. Higher order terms in the asymptotic expansion o//i,/2 can be simi- 
larly obtained, provided that V is sufficiently smooth. 
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Recalling p2|) . we see that for large the poles of ip can only come from the 
zeros of Wp. Substituting (fT4|) and (fTSj) into (fT3|) . we see that 

(19) Wp = (e-4^/=^/.r(p) + h,{p)] 
where 

(20) /ii(p) - (^>^y(i)~5'i(i)^ (^^^^n-i) + 32(i)) 

(21) /i2(p) =4ip3 + 2iVi(V^(l)-l^(-l))p^/' + 0(p2) 
and 

(22) /i3(p) = -2^^ + 0(1) 

Proposition 6. In the generic case when hi ^ 0, Wp has infinitely many zeros in 
the left half plane. Their asymptotic behavior is 

( ^i^k'^ -T,k log k + ayk + o{k) , V{1)V{-1)^0; 

(23) p= < -Vfc^- ^k\ogk + b^k + o{k), exactly one ofV{±l) is zero; 

[ ^iL^e -^-k\ogk + c,k + o{k), y(i) = v^(-i) = 0. 

where A: € N and k — )■ oo, and ay,by, Cy are constants. 



Proof. The equation Wp = reads 

<^*' - 

A simple analysis shows that this can only happen if p is near the negative imaginary 
line with p ^ —k^n^i/A where fc g N. We let p — ~i{ki:/2 + zY and rewrite ((24|l 
in terms of z: 

^=4z^"H" /^i(-W2 + .)^) , 

Recalling ((20)) and pT|) . we easily see that the right hand side of the above 
equation is contractive for large k. 

One can find the asymptotic behavior of z by iteration. First assume V {1)V {—1} ^ 
0. It is easy to sec that 

fe2(-»(W2)^) _ TT^k^ ^ 

~/ii(-i(fc7r/2)2) " 4F(1)F(-1)^ ^ ^ 
Therefore z ~ — ilogfc. Further iteration implies z = —ilogk + dy + o(l). 

Similarly, if exactly one of V{ztl) is zero, then z — — |ilogfc + by + o(l). If 
V{1) = V{-1) = then z = -|i log fc + c„ + o(l); Eq. ([23]) follows. 

□ 

The above analysis shows that all zeros of Wp for large p are in the left half 
plane. Thus we have 

Corollary 7. There are only finitely many bound states (this, of course can be 
simple shown by standard spectral techniques). 

We may now proceed to consider the order of these poles as well as their residues. 
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Proposition 8. The poles of "ijj for large p are simple, and the residues grow sub- 
exponentially. The residues of 1/Wp grow at most polynomially. (In fact, they 
grow exactly polynomially, since the asymptotic expansions in il9\) - [^) are differ- 
entiable.) 

Proof. Recalling we notice that 

Wp = w;{pk)ip~Pk){i + o{i)) 

where it can be easily checked that Wp{pk) 7^ 0. Then 1/Wp{pk) grows at most 
polynomially, and this together with the bounds on y± (see Lemma (HI) show that 
the residues of i> are bounded by o(elR<=v^l(kl+2)) □ 

The polynomial growth of residues, along with the analyticity of ip, show con- 
vergence of the sum in ^ as well as its Borel-summability. 

3.4. Asymptotics of We will show that V' has sufficient decay to allow for 
inverse Laplace transform as well as the desired bending of contour leading to 
Borel summation. First we rewrite (|10p as 

(26) - iWp4j{x,p) = y^{x) y+{s)ijo{s)ds - y+{x) y_ (s)V'o(s)ds 

Jm J-m 

assuming supp^/^o € [—M,M]. 

Lemma 9. y±=0 (^^e2|i^'=v^l(l^l+2)) for large p e C. 

Proof. We will prove the lemma for 2/+ using matching conditions. The proof for 
y^ follows analogously. 

The result is obviously true for x > 1, where y+{x) — e^^^*^^. 

For — 1 < x < 1, we have 

y+{x)^Cifi{x) + C2f2{x) 

where 

^i^^n^^ (/2(i)-/2(i)) 

C2^-V^—— (/l(l)-/{(l)) 

^f;p 

and Wf.,p = fi{x)f^{x) - f[{x)f2{x). 

It is easy to see, using Proposition [U that Ci,2 = 0(y/p). The bounds follow 
from dUl) and fTS]) . 

For X < —1 we have 

y+{x) = Cse^^^ + Cge"^^^ 

where 

C5 = ^e^(Vpy+(-i) + ^^v'+{-i)) 

C, = ^e-^(^j;+(-l) - Viy'+i-l)) 
The result is shown by estimating y^{—l) and ?/+(— 1) with Proposition [5l □ 
Lemma 10. p-^Wj,i){x,p) = O (^e2|R°v^l(l^l+*^+2) j for large peC. 
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Proof. A straightforward estimate from (|26|) and Lemma [9l 



□ 



Lemma 11. (i) There exists a set of curves Pk{s) parameterized by s € [0, 1] with 
p(0) on the negative imaginary axis, p{l) on the negative real axis, |pfe(s)| > k, so 
that 1/Wp is bounded uniformly in k by a polynomial in p along these curves. Here 
k € N can be chosen to be arbitrarily large. 

(ii)Moreover, 1/Wp is bounded by a polynomial in the region {p : argp G ct] 
and \p\ > Pa} where —n < a < —n/2 and Pa > depends only on a. 



Proof. We rewrite as 



W 

' ~ P'hsip) 



-ip 



We only need to show that 



on a chosen set of curves. 



h2{p) 



hi{p) 
h2{p) 



> 1 



Recalling the asymptotic expressions for /11.2, we have ^'iC' ^ '/4) 
co,n are constants. 



2_2 ■ 

/i2(fc^7r2j/4) 



cofc" where 



Let Pk{s) = — i(fc7r/2 — iargco)^(l — is^ and we have — 4-^— ipfe = {2ki:i 



jargco)(l — is) — (2/c7r — argco)s + ikiri 
Thus for s e [0,1/yfc] we have 



zargCQ. 



hi{p) 
h2{p) 



for all /c G N, while for s G [1/Vfc, 1] we have 



h2{p) 



2Vfe7r 



for aU fc > 0. 

The second part of the lemma follows from the above inequality since k may be 
taken to any large real number. Note also that Re(— -^Z— zp) = (k-K/2 — ^ argco)s. 
It is easy to see that pk (s) also satisfy the other conditions specified in the lemma. 

□ 



Collecting the above results we obtain 

Lemma 12. ijj{x,p) ~ O ^e2|ReV-ip|(|a;|+M+2)~j large p, in any given sector 

argp G [— TT, q;],|p| > Pa where —tt < a < —■n/2 as well as along curves Pk{s) as 
shown in the previous lemma. 

3.5. The inverse Laplace transform. To obtain the transseries of ?A from our i/), 
we take the inverse Laplace transform, lao-ix B''*4'{p)dp and push the contour 
into the left half plane. We will justify this procedure in this section. 
First we rewrite ([7]) as an integral equation 



y = T{Vy + iipo) 
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Figure 1. Curves Pk{s) passing between poles (plotted with 
square barrier potential) 



where 



(27) T{f)ix) —^e^-^P^ e'^^^^^ f{s)ds 



1 -''~^px / e^-'P'f{s)ds 



We further let y{x) — T{iipo){x) + p~^/'^h{x) and rewrite the integral equation 

as 

(28) h=p'"''T{y ■T{i4,^))+T{Vh) 

We start with a simple observation. 

Remark 7. e^^^^^ is bounded in the region := {p d C : —tt/2 < argp < 
tt} 1J{p G C : —Imp > const (Rep)^ . (In the following, we will choose const = 1/9) 



Note that 



1 



Rc{- y^-ip) = -— y/Rep + \p\ {Imp - Rep + \p\) 

Zimp 



We denote fi = sup^gQ |e ^ 

Lemma 13. Assume f and g are locally bounded functions and fg is compactly 
supported, with supp(fg) G [—b,b] where b> 0. Let a > b be an arbitrary number, 
fl' = fl \J{p e C : \p\ > py > 1}. We then have 

\Ti fg)\ < ^ ^P^g[-b,&] IsWI I 



y/Pv 

where \\f\\ suppgs^,^^g[_„ \f{x,p)\ 
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Proof. By ([27| we have 

(29) \T{fg){x,p)\^-^ [''\e-^-\\g{u + x)\\f{u + x)\du 



2|VPl 
1 





le'^'^Wgiu + x)\\f{u + x)\du 

-b 



^ ^' I f Miff ^ ^ ^^Pa:e[-b,b] 



J -b \/Pv 



□ 



Lemma 14. For compactly supported and twice dijferentiable ^po, we have 

T{iipo){x) = -'tpoix) + -^Gi{x,p) 
p p-^i^ 

where \Gi{x,p)\ < 2Af sup 1-00 I suPsg[Q^^,/_|_|3.|] |e^^~*P^|, assuming supp £ {—M,M). 
Proof. This is shown by repeated integration by parts to (|27|) . Note that 

(30) T{ii,o)ix) = -M^) - ^e^- r e-^'i;'o{s)ds 

- -^e-^^ r e^^^'ois)ds 

J-M 



1 



M 

e^-'P'ij'^{s)ds 



□ 

With the above lemmas, we have 

Proposition 15. Let Qq = ^[J{p £ C : \p\ > py} where py — 9(sup^g[_i ^(2;) + 
^ + 1)^. Let xi > be an arbitrary real number. The integral equation i28\) is 
contractive in the space of functions analytic in p € Vlo equipped with the sup norm 
11/11 = ^^Pp€no,xe[-~xuxi] \fix:P)\, within a ball of size 

2/x sup |F(x)| +2Msup|V'ol sup |e"^^"| 

a:e[-l,l] se[0,A/+|xi|] 

In particular, the solution h is bounded as xi ^ 00 ifKep > 0. 5*66 Fig. \3. 51 

Proof. The estimates of p'^^^T{V ■ T(i^o)) follow from lemma [H] and [T31 with 
f = V , g = ipo and g = Gi separately. 

The contractivity of T follows from lemma [T^ with f = V, g ^ h. Note that 
analyticity in p is preserved by T and convergence in the sup norm. □ 

We therefore have the following results. 

Proposition 16. (i) ip (as in section [3. 1\} has the following decomposition: 
i'{x,p) = -tpoix) + -^G2{x,p) 
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where G2{x,p) is hounded in p ^ ^o, x G [—xi,xi] where xi > is arbitrary. 

^(x, t) = Vo(x) + r^^^ ^^^e^Up 
is the solution to ([7p. Here oq > is a constant. 

Proof. We only need to show that in the solution tp is identical to the solution y 
obtained in this section, the decomposition of which has already been shown. Part 
(ii) then follows immediately from properties of the inverse Laplace transform. 
To this end, note that the general solution to (O can be written in the form of 

ygen{x,p) = ■tpix.p) + ci{p)y+{x,p) + C2{p)y^{x,p) 

where y+ and y- are the homogeneous solutions defined in section lSTT] (with a slight 
abuse of notation). This implies 

y{x,p) = i^{x,p) + ci{p)y+{x,p) + C2{p)y-{x,p) 

where y{x,p) is the solution obtained earlier in this section. Since in the region 
{Re{p) > Pv, X < 1}, y+ is unbounded and y_ is bounded, and in {Re{p) > Pv, 
a; > 1} y+ is bounded and ?/„ is unbounded, while both y and -0 are bounded 
(the boundedness of ■0 follows easily from (|26p ). we must have ci = C2 = in 
Re(p) > py. Thus and y coincide in Re(p) > p^, and also in f^o by uniqueness of 
analytic continuation. 

□ 

Lemma 17. In the expression 

j_ n-- b^,.up 

2-Ki / r)3/2 P 

^T^'' Jao-loo P 

we may deform the contour to one which goes from —oo below the real axis, turns 
counterclockwise around the origin and goes towards —oo above the real axis. In 
the process we collect all residues from all the poles in the left half plane. 

Proof. The deformation of the upper half of the contour is justified by Proposition 
1161 since f2 contains the second quadrant. 

In the third quadrant, recall that Pk{s) = —i{kTT/2 — ■|argco)^(l — is)^ and 
Re(v'— ipfc) = {kn/2— iargco)s. Thus Repfc(s) ^ const. k^s and Ke{y/—ipk) = 
0(Repk{s)/k) for all s > 0. Therefore we may choose part of the curve pfc(s) where 
s S [l/fc,l] and join it with a curve in {p S C : —Imp > (Rep)^ /9} f]fto, say a 
vertical line downward to infinity. These two curves, along with the one from below 
the real axis to the origin and lower half of the original contour, surround all poles 
in the third quadrant as fc — > oo. Decay along the Pfc(s) curve is ensured by the 
term e^*, since eP'=(^)*±V-«P'=(^)A^o = 0(6"'=*) for arbitrarily large Mq. Note also 
that the length of Pk{s) is of order k^. □ 

To prove Theorem [I] we further write ipo{x) — J^, ^ipoix) and combine it 

^^^^ 2^ Ic ^^3/2^'* eP*^dp, where the contour of integration is the horizontal part 
around the negative real axis described above. This contour can be deformed to 
"0 to — oo" in an upper and a lower sheets of the Ricmann surface, which yields a 
Borel-summable power series in t^^/^. 
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Figure 2 . A sketch of the contour deformation used. The shaded 
area is the contractivity region. 



3.6. Connection with Gamow Vectors. Classically, Gamow vectors are ob- 
tained as solutions to ([8]) with "purely outgoing boundary conditions" as x ^ ±oo. 
In our case, this means such a solution (after rescaling) equals y+{x) for x > 1 
and a nonzero constant times y~{x) for x < 1, y± being as in section (|3.ip . The 
existence of such a solution, therefore, is equivalent to the linear dependence of t/+ 
and ?/_ (cf . Lemma , which in turn is equivalent to the vanishing of the Wron- 
skian: Wp = 0. Thus the 7^ found from the poles of tp are exactly the resonances 
corresponding to the Gamow vectors, a constant multiple of y+. The latter are is 
easily seen to be multiples of the residues of for example by simplifying pB]) : 



-iWpipix.p) ^ cy+{x) I y+{s)%l)n{s)ds ~ y+{x) I cy+{s)%l)n{s)ds 

JM J-M 

3.7. Proof of Proposition [2j 

Proof. This follows straightforwardly from Lemmas [TH] and 15.11 after extracting a 
suitable number of poles from ip. All poles are simple, and the contribution of a 
pole of residue and position pk is 

Jk\t) = Tk / —p 

J a Vp - Pk 

The representation of Jk in terms of special functions is perhaps most conveniently 
shown by solving the first order differential equation it satisfies, and determining 
the free constant from the asymptotic behavior m Pk- D 
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4. Example: the square barrier 

Here we take as a simple example the Schrodinger equation with a square bump 
potential V{x) — X[-i,i]j X being the indicator function. (One of few cases where 
explicit solutions exist.) 

Aie^^^ + Aae-^^^, x ^ -1; 
= <( Agev^i^^ + ^46-^^^^^, -1< a; < 1; 
^-V^x^ x > 1. 

y^{x)^{ Bie^-^i^^ + Bae-^-"!^^, -1< x < 1; 

where the coefhcients Aj, Bj are determined by matching solutions at the endpoints, 
±1. For example, 



The other coefficients are similar (and obtained in a similar way) and we omit them. 

It follows that the Wronskian Wp has an explicit expression 
(31) 



{i + 2p - 2^\/ i+p) — i ~2p— 2^\J i + p 



2^/TTp V 

We may find the asymptotic positions of the resonances by iterating 



.„„x 1 , f i + 2pk + 2VPfe + Pk 

(32) ^fc = 1- log 



4i \i + 2pk-2Jpk^/r+p^ 



where pk = — i(fc7r/2 + z^)^. 

We also calculate the residues of 1/Wp by differentiating (PT|) : 



, ^/Pk{^ +Pk){i + 2pk - 2^/pl^i+pk) 1 

i/ Wp ~ 



Here we calculate the positions and residues of a series of poles using the above 
formulas and compare them to the asymptotic behavior — 7rfclog(7rfc) — i7r^A:^/4, as 
in Proposition [6l Then we plot these poles together with a density graph. 



The asymptotic pole location formula gives (increasingly) good accuracy starting 
with the 15th pole or so, where it predicts the position —181 — 555i, whereas the 
exact value is about —180 — 532z. 

The first resonance, the one closest to the imaginary line (in p-plane), may have 
a visible effect on the wave function ■(/' even if this resonance does not correspond 
to a (long-lived) metastable state. We will demonstrate this phenomenon, as well 
as the computational effectiveness of the Borel summation approach, using (near-) 
optimal truncation, see i)5.11 on the example of the square barrier potential, where 
we choose the initial condition to be V'o(a;) = x[~5; 5] simplicity. 

In our example, the first pole of 1/Wp is located at po = —1.70018 — 0.805871i. 
This can be found by standard iterative arguments. 
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-280 -260 -2W -220 -200 -180 -160 -IW 



Figure 3. Density graph of 1/Wp. Dark dots indicate poles cal- 
culated from the asymptotic formula. 



We will demonstrate the effect of this pole in the region x > 1, where (cf. (fTO|) 
ip{x,p) = — / y-{s;p)ds 

and 



Wp 



(33) tp{x,t) = -—l — / y^{s;p)dsdp 

= loo 7-. ^-(^^^^'^^'^^ - 2^ X 1. y-is;p)dsdp 

+ hm ((p-po)^(a:,p))e«*(l + o(l)) 

for large t. 

We may calculate the power series by expanding ip{x,p) near p — and using 
Watson's Lemma. For instance, for a; = 8 we obtain the series 

(34) 

(0.735266 + 0.735266i)47:r - (12.3883 - 12.3883z)47t - (98.5277 + 98.5277i)47^ 
+ (471.935 - 471.935i)-J- + (1429.08 + 1429.08i)-i- 
- (2690.72 - 2690.72i)^Y^ - (4000.95 + 4000.95^)^^^ + OilT^) 

Note 8. Taking x = 7, the contribution of the first resonance to the power series is 
visible, about 3%, for t — 7, and evidently decreases exponentially thereafter. The 
overall precision increases rapidly with t, if x is fixed or does not increase faster 
than t, unlike most direct numerical calculations. 
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5. ECALLE-BOREL SUMMATION, EXPONENTIAL ASYMPTOTICS 

An expansion is Borel summable if, by definition, it is the asymptotic expansion 
of the Laplace transform £ of a function which is real-analytic on R"*", exponentially 
bounded, and which has a convergent series at zero, in (ramified) powers of the 
variable and possibly logs (Frobenius series). The Borel summation operator, CB^ 
is essentially CSC^^ where applied to a series is understood in the formal 
sense, as the term-by-term transform, and S is convergent summation. Since CB is 
conjugated with usual summation, which commutes with virtually all operations, 
the same is true for CB. Uniqueness of the Borel sum stems from uniqueness of 
the sum of a convergent series. Borel summation is a canonical extension of usual 
summation. The decomposition of a function in a Borel summed part and a sum of 
exponentials, when possible, is also unique and canonical, see ^and for a detailed 
analysis e.g. [5J. 

5.1. Borel summation and least term truncation. Borel summation allows 
for exponentially accurate calculations of the associated function by truncating the 
series near its least term. We first briefly explain the reason and refer to [7] for 
more details. The accuracy, even for t not so large, is illustrated in 21 

Here we consider a Borel summed scries of the type that intervenes in our prob- 
lem, namely let 

/•oo 

(35) f{t) = / e-P'F{p)dp 

Jo 

where F{p) = g{p^);Ile(3 > and g is analytic at the origin and mcromorphic in 
C. The asymptotic behavior of / is, by Watson's lemma, of the form 

oo 

(36) /(i)^/(i)=^c,t-^-^ 

k=0 

We let Sn{t; f) be the truncate of the power series, up to the power t": 

(37) Snit-J)^ ^fe^""" 

0<feRc,3<n 

We want to estimate the error by calculating the function from its power series 
by optimal truncation, or truncation to the least term. For a series in which the 
coefficient of t'' grows roughly like b\/c^, means using the truncation Sct{t; /). For 
example using Stirling's formula we see that for F{p) = l/{p — po), the general 
term of / grows like fc!/|po|'^ and its least term is near n — t|poh this location 
of the least term is the same regardless of the nature of the singularity, for all 
algebraic-logarithmic type of singularities. 

We show that the error in approximating / by St{t; f) in this way is of the same 
order of magnitude as this least term, which is exponentially small in t. In similar 
contexts this is known quite generally for Borel summed series, see [7] and references 
therein. Least term truncation provides a practical way to calculate functions with 
very high accuracy even for t of moderate size. In spite of the generality considered 
in [7], our case is not covered (because of ramification at zero). Instead of describing 
the adaptation of that proof, for convenience of the reader, we provide a complete 
argument in our case. 
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Lemma 18. Assume (3 G (0,1) and H is analytic in disk Byi2- Let Hm be the 
maximum of H on a disk of radius A £ (^1,^2)- Then 

41/0 



(38) H{p^^)expi~tp)dp^S^^+^^p{t,H{p'')) = 0^ ^n+i^^^p+p+i 
Proof. Indeed, by using Taylor series with Cauchy integral remainder we have 

(39) / i7(/)exp(-ip)dp - SniHip^)) 
Jo 

where C is a circle of radius A. If is the maximum of H on C, We have 

and the result follows. □ 

Lemma 19. Assume j3 £ (0, 1) and let F{p) = [p^ ^ Po)^^ where po ^ IR+ and let 
f = CF. Then, 

(41) i/w-5,.,.,,(u-)i< j.';:^:,f.;j^^, 

Proof. Writing 



(42) 



1 " pjP 



P-Po Up'o PoiP^-Po) 



this follows from straightforward integration and estimates. □ 

5.2. A class of level one transseries. We only need an especially simple subclass 
of transseries, exponential power series of the type 

00 

(43) />)=E'^"^'**"'/fcW 

where fk{t) are formal power (integer or noninteger) series in where, for dis- 
ambiguation purposes, the real part of the leading power of 1/t in fk{t) is chosen 
to be 1. Agreeing that no fkit) is exactly zero and the are distinct, it is required 
that the exponentials e~'^''* are well ordered, in the sense that Re(7fe) > Re(7fe') if 
k > k' , and every Re(7fc) has a predecessor, the smallest Re(7j) greater than it. In 
our context the sets {j : Re(7j) = Re(7fc)} turn out to be finite. 

The transseries / is Ecalle-Borel summable if (a) fk{t), k £ N are simultaneously 
Ecalle-Borel summable (in fact, simply Borel summable, in our case), and (b) upon 
replacing each fk{t) by its sum, the resulting function series is uniformly convergent. 
We give precise definitions in ^ Transseries and Ecalle-Borel summability were 
introduced by Ecalle in the 1980s and there has been substantial development since. 
For an elementary introduction see [5] . 

The transseries is (Ecalle-Borel) summable if for some T > we have the fol- 
lowing. 
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(i) fk{t) are simultaneously Borel summable, that is there exists a T so that 
fk{t) are the asymptotic expansions for large t of Laplace transforms, 

/■oo 

(44) fk{t)^ Fk{p)e-P*dp=:CBMt) 

Jo 

where 

(ii) Fk are ramified-analytic at zero, and real analytic on M+ with the uniform 
bound \\Fk{p)\\ < CkcM^. 

(iii) For some e M we have have \Fk{t)\ < Cke^^P^ . 

(iv) The series 

oo 

(45) El^^'^'^l^'^ 

k=Q 

converges for some T > (and thus for all t > T). We recall that, by convention, 
Fk{p) = Cfc(l + o(l)) as p ^ 0, where Ck 7^ 0. 
Therefore, the sum 

00 00 00 

(46) f = CBY,t"''e~'"''fk{t) :-^e-^'=*r''/:S/fc(i) = ^e-'^^*t"'=/fe(t) 

A:=0 fc=0 fc=0 

converges absolutely for t > T. 

The operator CB is a proper extension of the Borel summation operator. In par- 
ticular, it allows for non- accumulating singularities on the axis of summation, in 
which case analytic continuation is replaced by Ecalle's universal averaging. Super- 
exponential growth of F of a controlled type is allowed, using Ecalle's acceleration 
operators. 

With these extensions, Borel summation is an extended isomorphism between 
series, or more generally transseries, and a class of functions {analyzable functions), 
commuting essentially with all operations with which analytic continuation does. In 
this sense, Ecalle-Borel summable transseries substitute successfully for convergent 
expansions; in particular the Ecalle-Borel sum of a formal solution of a problem 
(within certain known classes of problems such as ODEs and PDEs) is an actual 
solution of the same problem. It is known that the fundamental decaying solution 
of a nonlinear differential equation at a generic singularity is given, uniquely, by 
Borel summable transseries [3]. 

5.3. Uniqueness of the transseries representation. In the same way as the 
asymptotic power series of a function, when a series exists, is unique one function 
can only have one transseries representation, if at all. We sketch a proof that a 
representation of the form of a given / is unique. We assume of course that 
the transseries are in canonical form, as explained above. By linearity, it suffices 
to show that if / given in (|46|) is identically zero, then all fk, and thus all fk are 
identically zero. We assume by contradiction that some fk are nonzero. Since the 
Re(7/c) are well ordered, cf. ^5.2i we choose the largest Re(7fc) such that fk ^ 0. 
There are only finitely many 7^ with the same Re(7fc), cf. again ^5.21 We can 
assume without loss of generality that these As have indices 0, ...,n, and assume 
that we have ordered the terms in the transseries so that Re7i < Re7i+i for all i. 
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We write 

n oo 

(47) / = ^e-^'=*r'-/fc(i)+ E e-'"'H^^ fk{t) 

k=0 fc=n+l 

Note that for any e > small enough we have 



(48) 



J2 e-'^^^+^k'-' fkit) 



k=n+l 



< const 



as T — > CXD 



since Re(7o) < Re(7„+i). Dividing (gT]) by e"^"'^''* we get 



(49) 



k=0 



'r'fk{t)^o{i), (t^oo) 



where ak = Ini7fc. For each k we choose (3k to be the smallest power of p (in 
absolute value) with nonzero coefficient, Ck in the expansion of Fk- Of course, if 
all coefficients in the Puiseux series of Fk vanish, then Fk vanishes near zero, and 
thus everywhere by analyticity. We arrange that there is no k such that Fk = 0. 
Then, by Watson's lemma, Fk = CkT{Pk + l)i~''''~^(H-o(l)) for large t. We choose 
the largest (3j, in the sense above, and divide by r(/3j + l)t^^'^^^^^ . We get, by 
Watson's Lemma, 



(50) 



RcfJj=Rc[3k;k<n 



o(l) 



where 9k = Im/3/c. We now prove a lemma in more generality than needed here, in 
view of future generalizations to time dependent potentials. 

Lemma 20. Assume X^fc^o I'^fcP °° '^"'^ 



(51) 



/(t)=^Cfee-*"-4-^''- =o(l) 



fc=0 



where ak, 0k G K. as t ^ oo. Then f{t) = 0. 

Proof. We first look at the simpler case where all 9k — 0; as we shall see, the general 
case is similar. We see, by explicit integration and dominated convergence, that for 
large to and t — + oo we get from ([54| that 



/ |/(.)pd. = 5]kpt + o (y 

•^*o k=0 \k=0 



Ck\ 



o{t) 



(52) 

which is only possible if 

oo 

(53) Y\^k\' = 

k=0 

To generalize to the case 9k 7^ 0, we simply note that ([54| implies 

00 

(54) /(e^) = Y Cke"""-'" e-^'"' = o(l) 

fe=0 
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OO we have (e.g. by integration by parts) that, for 

^e-'""'e-'''''(l + o(l)) 


□ 

Borel summation and usual summation: the underlying isomorphism. Further- 
more, there is the following important point. When a Borel summable transseries 
of a function exists, functions and their transseries have the same properties. That 
is, there exists an extended isomorphism between transseriable functions and trans- 
series similar in many ways to the one between germs of analytic functions, and 
their local convergent Taylor series regarded as formal algebraic objects. This latter 
isomorphism is so flawless that we do not distinguish notationally a convergent sum 
as a formal sum, from its sum as a function. Borel summation is a proper extension 
of usual summation, carrying further these isomorphism features. 

The isomorphism, provided by Ecalle-Borel summability, which recovers the 
function from its transseries, justifies the usage of the term complete asymptotics. 
Borel summation is a canonical way to sum factorially divergent series, cf. also [5]. 

Independence of method. Finally, the nontrivial terms in the transseries of a 
function can be exhibited by many other exponential asymptotic techniques some 
of which having of substantial calculational value, such as hyperasymptotics, a set 
of methods improving and refining optimal truncation of series, cf. [2], [5], [17j . 
and references therein. 

In the language of generalized Borel summability, the wave function asymptotics 
is given in all amplitude regimes by an Ecalle-Borel summable transseries, valid for 
t > 0, and this transseries turns out to rest on a Gamow vector decomposition. 

Note 9. Sometimes a given series can be Borel summed with respect to different 
powers, or more generally functions, of the variable. For instance, 

(56) r e-P'e-'P'/^p^ /"e-P*'— ^ 

^ ' Jo Jo 2^{i-p) 

Since the integrals are equal, they have the same asymptotic series for large t; both 
integrals are Borel sums of the same asymptotic series, on the left interpreted as a 
series in 1/t, while on the left it is thought of as a series in \/t^ . 

Using the connection with Gevrey asymptotics, [S], it is easy to see that a series 
has a unique Borel sum, with respect to any variable in which it is Borel summable, 
even when allowing for ramified-analytic functions. Ramified-analytic functions 
are real analytic, and near p = of the form F{pi, ...,pm) where pj — p""^ \ogp^^ , 
Re(ai) > 0, and F is analytic at . It is an easy exercise to show that e~* , Rea > 
cannot be represented as a Laplace transform of a ramified analytic function. 

But beyond ramified analyticity uniqueness of the representation as a "contin- 
uum" (to use physics terminology) plus exponentials docs not hold. We have, e.g., 

/•oo 

(57) e-* = 7r-i/2 / p-^/^e'^/^e-P'" dp 

Jo 

which is a continuum type integral. 

The Borel sum gives consistent results, and remaining exponential terms are 
uniquely defined. 



as s — + oo and that, still as s 
9 ^ 0, we have 




e 
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Ramified analyticity of ^jj follows from the formulas in 1Q\ and However, 
the techniques [8] are more involved and, along those lines, there appear to be 
significant gaps in estimates leading to a mathematical proof of Borel summability 
(which, in fact, is not the intention of those works). The purpose of the analysis 
in [8], [16] and related literature is different: the extension of a spectral- like theory 
and a "spectral calculus" beyond the continuous spectrum. 

5.4. Analytic potentials. Since the wave function is the solution of a PDE which 
mixes space and time information, finding the detailed time behavior of ij: is contin- 
gent on detailed information about V{x) and '00(2;). It seems likely that generalized 
(multi-) summability of the large time (trans-)series of '0 should hold whenever V 
has a multisummable transseries as well. In this paper though we consider poten- 
tials analytic at infinity and with sufficient decay. For simplicity, we write —ip ~ e^. 
Equation ([8]) reads: 

(58) y"{x) = {V{x) + e^) y{x) 

From the form of the Green's function, it is clear that the analytic properties of 
the Green's function at p = follow from those of the Jost functions (defined as 
in (O: is the solution that behaves like e~^~*P^ as a; oo, when p G IR+). 
We analyze potentials of the form V{x) = a/x"^,m S N, and we discuss how 
essentially the same arguments would work for any potential which is analytic at 
infinity and ©(x"™). The value m = 2 marks in a sense a threshold, making the 
transition between convergent and divergent expansions in energy at the bottom 
of the continuous spectrum. For m ~ 1,2 the equation can be solved in terms of 
simple special functions; the slow decay in the case m — 1 implies that zero is an 
accumulation point of poles; no convergent Frobenius expansion is possible. 

Proposition 21. For m > 2 the Jost functions have convergent Frobenius expan- 
sions in e (series in ramified powers of e and eloge and are bounded by conste^'^^^''^ 
uniformly in a sector arge € (— 37r/4, 37r/4). 

For V{x) — a/x^, (l58|) is solved by Bessel functions; the solution that decays like 
e~'^^ as a: ^ oo is given in terms of the Bessel function K as 

?/+ = y/2ex/TrKa{ex); a = \J a +1/4 

For small e and fixed a;, y"*" has the form 

(59) Ci(a;)ei+"^i(e) + C^ixy-'^ A^ie) 

with A\ , A2 analytic. For m > 3, there are no ramified powers of e in the expansions, 
but all powers of eloge intervene. 

Proposition 22. For fixed x and m > 3, the function s{x; e) is of the form 
G(e, eloge) where G{u,v) is analytic for small (it, w). The Jost functions are 
bounded by conste^'^^^^^ uniformly in a sector arge S (— 37r/4, 37r/4). 

Proof. The question is the dependence of the Jost function in e, for small e e C. 
As a mathematical question, this is a connection problem: the definition of the 
Jost function is given in terms of the asymptotic behavior as a; — > 00 while the 
analyticity properties in e are sought globally in x. 

It is convenient to treat this problem by Borel summation once again, this time 
in X, to transform it into a pure analyticity question. We analyze the Jost function 
given, for e > 0, by y{x) ^ e^'^^(l -I- s(a;;e)) where s(a;; e) is an o(l) power series 
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in l/x, as a; — > oo. It is easy to see that such a solution (whose existence is known 
and also follows from the argument below) is unique. 

For Borel summability, we have to extract s. We thus write y{x; e) = e~'^^{l + 
s{x; e)) and obtain 

(60) s" - Its - V{x)s = V{x) 

To simplify even further the presentation we take to = 3, but there is nothing 
special about this choice, and the extension to other values of to is immediate. 

We inverse Laplace transform ([50]) (the legitimacy of which is justified "back- 
wards" by showing that the Laplace transform of the solution of the thus obtained 
equation solves (|60[) . which has a unique small solution) and obtain 



(61) H^^Av'At^' "^'^ 



2^ 2 qiq + 2e) 

where VF is the definite antiderivative of F which is zero at zero. With the change 
of variable r = eg, H{q) = F{t), we obtain 



3 Fir) 



— +^^.(r + 2) 

We look for a solution which are 0{e^T^) for small t. Consider the space B of 
functions of the form /(t) = t'^G{t) where G is analytic for, say, |r| < 1 with the 
norm ||/|| = sup|T-|<i |G'(t)|. We see that this is a Banach space, and eq. ([6T|) is 
contractive in B. It is also unique in the space of functions of the form t'^G{t) with 
G defined in £^[0, 1]. The solution of ((6T|) is unique, and analytic for small r. As 
a differential equation this reads 

(63) F'" = 



t{t + 2) 

The argument above, or Frobenius theory, shows that (I63|) also has a unique solution 
which is of the form ie^T^(l + o(l)) for small r. The solution is obviously analytic 
for Rer > —2, since there the only singularity of the equation is r = 0. 

By standard ODE asymptotic results (l^ we see that any solution of ((63)) is 
uniformly bounded in C by 



(64) G\ 



2/3^3|r|i/3 



r ' e 



for some C. This ensures the necessary (sub)exponential bounds for taking the 
Laplace transform. 

On the other hand, we look for solutions of (l63|) in the form 

(65) F = e^F2 + J2^'PAr) 
The functions Fj satisfy the recurrence 

(66) f;'U = ^^ j>3 

■' T[T + 2) 

With our initial condition, we get F2 = t^/2 and 
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For now, we take r in the right half plane, H. It can be checked by induction that 
Fj are analytic in H and at zero, and 

(68) l^.l<4^ 

It follows that the series converges uniformly on any compact set in H. More- 
over, we see that the function series 

(69) H{q)^'^+Y,^=F,{q/e) 

also converges uniformly in q on any compact set in H. The Laplace transform of 
H reads 



poo 1 poo 1 

/ e-'^^H{q)dq = _ + ^ e""- F,{q/ e)dq Mx; e) 

"'0 X Jq X 



(70) 

j>3 ~ j>3 

where, once more, the interchange of summation and integration is justified by the 
bound (I68p . We fix a;, drop it from the notations, and note that in the last sum 
we have \fj\ < const{j\)~^ . We thus need to study the analyticity of fj. We claim 
that fj{e) — Gj(e, eloge) where Gj{u,v) is analytic in small (w, w). Dominated 
convergence ensures that the integral on the left side of ((70)) is of the same form. 

We will use the following Lemma which applies at the other end of Watson's 
Lemma setting. 

Lemma 23. Assume F is bounded on 0,M and analytic in a sector {z : \z\ > 
i?; arg(z) e {o.,b)} and f is of the form z^G{z~^, log z) where G{u,v) is analytic 
in the polydisk {{u,v) : \u\ < Mj~^, |u| < Mj~^}, let Mi > M and consider 



oo 



fis) = \ e-^^F(j>)dv 
Jo 

Then 

(71) f{s)^s-''-^H{s, slogs) 

where H(u,v) is analytic for small u,v. (In a very similar way, the lemma could 
accommodate for fractional powers of z~^.) 

Note that the convergence of the series in 1/x, x~^ log a; entails that F extends 
analytically on the Riemann surface of the log in a neighborhood of infinity. 

Proof. The proof is elementary. Let M > Mi, and first note that J^^^ e~'^PF{s)ds 
is entire, and we only need to consider the integral from M to infinity. Since for 
some C > and all {l,j) > (0,0) we have 

POO 1 

(72) / e-">^-"^' log^ {p)dp < constM^-"-^ log^ {N) - 
Jm s 

the series 

(73) F(p)=p^^Cfc,x-'=(x-ilogx)' 

k,l 

can be integrated term by term and uniform convergence easily entails that it is 
enough to show the property for a single term, of the form 

/"OO 

(74) Q(s)= / e'>^-"-'log^(p)dp 

JM 
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A finite number of integrations by parts, multiplications by s and differentiations 
in s brings ([73]) to 

(75) / e-'Pdp = e-^'ys =-+ entire{s) 

Jm s 

Undoing the operations above on tfie last expression in ([75)1 easily completes the 
proof. □ 

For large enough a we now write 
(76) 

poo pM poo 

Jo Jo Jm 

where we choose M large enough. 

The first integral is manifestly entire in e. For the second term, we have the 
following. 

Lemma 24. Fj{s) — s^Wj{s~^,s~^\ogs) where Wj{u,v) is analytic for small u, v. 

Proof. Induction from (j67p : The right side operations on the right side consist in 
multiplication by t~^{t + 1)^^, and three definite antiderivatives (from zero). It is 
sufficient to show that each of these operations preserves the structure above, aside 
from the leading order behavior which follows from straightforward power counting. 
Multiplication by r~^(r + 1)"^ clearly preserves the structure mentioned. 

(77) 



/ ^0 JM 

where the first integral is a mere constant, and M is chosen so that Wj is analytic 
for \u\ < 1/M and \v\ < 1/M. We then write Fj = Y.k<j+i.i '^j-M^'^v^ + Fji 
with u — s^^, V ~ s^^ logs and where we see that Fji = 0{s~'^ logs''). The sum 
contains finitely many terms, and for it the structure follows by explicit integration. 
For the second we write 

(78) 

iM Jm 

where the first integral on the right is a constant. For the second one, the structure 
for large t follows from term by term integration and straightforward estimates. □ 

The rest of the proof follows from Lemma [221 noting that, for fixed I, there are 
only finitely many terms in the expansion at infinity of / for which the total power 
of T exceeds —I. 

It is clear that all the arguments above go through if x^^ is replaced by x~"\m > 
3, except for ([M)) where the exponent will be p^/™ and the power of the prefactor 
changes. The bounds for the Jost functions follow immediately from the Laplace 
representation of s{x) and contour deformation. 

If the potential is analytic and 0(a;~™) at infinity, then the fimction H will, in 
general, have exponential order one, rather than fractional, and the bounds (|68p 
are "worse" , the power of the factorial becoming one. This can be shown similarly, 
using a roughly similar recurrence. In the 0(x^^) case, one would get recurrence 
of the form 
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where aj grow at most geometrically. The rest of the proof is roughly the same, 
but the details are more cumbersome. □ 
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